TREES, PARKING FUNCTIONS, SYZYGIES, AND 
DEFORMATIONS OF MONOMIAL IDEALS 



ALEXANDER POSTNIKOV AND BORIS SHAPIRO 



Abstract. For a graph G, we construct two algebras, whose dimensions are 
both equal to the number of spanning trees of G. One of these algebras is 
the quotient of the polynomial ring modulo certain monomial ideal, while the 
other is the quotient of the polynomial ring modulo certain powers of linear 
forms. We describe the set of monomials that forms a linear basis in each of 
these two algebras. The basis elements correspond to G-parking functions that 
naturally came up in the abelian sandpile model. These ideals are instances 
of the general class of monotone monomial ideals and their deformations. We 
show that the Hilbert series of a monotone monomial ideal is always bounded 
by the Hilbert series of its deformation. Then we define an even more general 
class of monomial ideals associated with posets and construct free resolutions 
for these ideals. In some cases these resolutions coincide with Scarf resolutions. 
We prove several formulas for Hilbert series of monotone monomial ideals 
and investigate when they are equal to Hilbert series of deformations. In the 
appendix we discuss the abelian sandpile model. 



1. Introduction 

The famous formula of Cayley says that the number of trees on n + 1 labelled 
vertices equals (n + l)™ -1 . Remarkably, this number has several other interesting 
combinatorial interpretations. For example, it is equal to the number of parking 
functions of size n. 

In this paper we present two algebras A n and B n of dimension (n + l) n_1 . The 
algebra A n is a quotient of the polynomial ring modulo a monomial ideal; and the 
algebra B n is a quotient of the polynomial ring modulo some powers of linear forms. 
It is immediate that the set of monomials x b , where b is a parking function, is the 
standard monomial basis of the algebra A n . On the other hand, the same set of 
monomials forms a basis of the algebra B n , which is a non-trivial result. 

More generally, for any graph G, we define two algebras Ac and Bq and describe 
their monomial bases. The basis elements correspond to G-parking functions. These 
functions extend the usual parking functions and are related to the abelian sandpile 
model; their number equals the number of spanning trees of the graph G. This 
implies that dim Ag = dim Bq is also the number of spanning trees of G. 

All these pairs of algebras are instances of the general class of algebras given 
by monotone monomial ideals and their deformations. For such an algebra A and 
its deformation B, we show that dim .4 > dim£> and the Hilbert series of B is 
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tcrmwise bounded by the Hilbert series of A. There is a natural correspondence 
between polynomial generators of the ideal for B and monomial generators of the 
ideal for A. However, these monomials are not the leading terms of the polynomial 
generators for any term order, because they are usually located at the center of 
the Newton polytope of the corresponding polynomial generators. The standard 
Grobner bases technique cannot be applied to this class of algebras. 

We also investigate the class of order monomial ideals that extends monotone 
monomial ideals. These ideals are associated with posets whose elements are marked 
by monomials. We construct a free resolution for such an ideal as the cellular 
resolution corresponding to the order complex of the poset. This resolution is 
minimal if the ideal satisfies some generosity condition. In this case, the numbers 
of increasing fc-chains in the poset are exactly the Betti numbers of the ideal. This 
resolution often coincides with the Scarf resolution. 

We discuss some results of our previous works on the algebra generated by the 
curvature forms on the generalized flag manifold. This algebra extends the coho- 
mology ring of the generalized flag manifold. For type the dimension this 

algebra equals the number of forests on n + 1 vertices. The algebras generated 
by the curvature forms are analogous to the algebras that we study in the present 
paper. This attempt to lift Schubert calculus on the level of differential forms was 
our original motivation. 

The general outline of the paper follows. In Section we define G-parking 
functions for a digraph G. We formulate Theorem 12 . 1 1 that says that the number 
of such functions equals the number of oriented spanning trees of G. Then we 
construct the algebra Aa as the quotient of the polynomial ring modulo certain 
monomial ideal. Elements of the standard monomial basis of Ag correspond to 
G-parking functions. In Section |3 we construct the algebra Bq as the quotient of 
the polynomial ring modulo the ideal generated by power of certain linear forms. 
Then we formulate Theorem 13.11 that implies that the algebras Ag and Bq have 
the same Hilbert series. In Section 0] we give two examples of these results. For 
the complete graph G = K n+ \ we recover the usual parking functions and the 
algebras A n and B n of dimension (n + l) n_1 . For a slightly more general class 
of graphs we obtain two algebras of dimension I (I + kn) n ^ 1 . Sectional is devoted 
to description of monotone monomial ideals and their deformations. We formulate 
Theorem 15.21 which implies the inequality for the Hilbert series. In Section [B] 
we describe a more general class of monomial ideals associated with posets and 
construct free resolutions for these ideals. Components of the resolution for such 
an ideal correspond to strictly increasing chains in the poset. In Section [3 we give 
several examples of minimal free resolutions. In Section|S|we prove general formulas 
for the Hilbert series and dimension of the algebra given by a monotone monomial 
ideal. Then we deduce Theorem 12.11 In Section we construct the algebra Cq 
and prove Theorem 19.11 that claims that the dimension of this algebra equals the 
number of spanning trees. Actually, we will later see that Cq is isomorphic to 
the algebra Bq- In Section ^]] we prove Theorem 15.21 Then we finish the proof 
of Theorem 13.11 which goes as follows. By Theorem 15.21 and construction of Cq 
we know that Hilb Ag > Hilb Bq > Hilb Cq termwise. On the other hand, by 
Theorems 12.11 and 19.11 dimple = dim Cq is the number of spanning trees of G. 
Thus the Hilbert series of these three algebras coincide. In Section ^2 we discuss 
some results of our previous works and compare them with results of this paper. 
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We mention a certain algebra, whose dimension equals the number of forests on 
ri + 1 vertices. This algebra originally appeared in the attempt to lift Schubert 
calculus of the flag manifold on the level of differential forms. In Section ^| we 
discuss a special class of monotone monomial ideals and their deformations. We 
give a minimal free resolution and subtraction-free formula for the Hilbert series 
of the algebra A and list several cases when it is equal to the Hilbert series of B. 
The appendix is devoted to the abelian sandpile model and its links with G-parking 
functions. 
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pleted, and by NSF grant DMS-0201494. The second author is sincerely grateful 
to the Max-Planck Institut fur Mathematik in Bonn for the financial support and 
the excellent research atmosphere during his visit in 2000. 

2. G-PARKING FUNCTIONS 

A parking function of size n is a sequence b = {pi, . . . , b n ) of non-negative integers 
such that its increasing rearrangement c\ < ■ • • < c„ satisfies c, < i. Equivalently, 
we can formulate this condition as | bi < r} > r, for r = 1, . . . ,n. The parking 
functions of size n are known to be in bijective correspondence with trees on n + 1 
labelled vertices, see Kreweras |Krew] . Thus, according to Cayley's formula for the 
number of labelled trees, the total number of parking functions of size n equals 
(n + 1)" _1 . In this section we extend this statement to a more general class of 
functions. 

A graph is given by specifying its set of vertices, set of edges, and a function that 
associates to each edge an unordered pair of vertices. A directed graph, or digraph, 
is given by specifying its set of vertices, set of edges, and a function that associates 
to each edge an ordered pair of vertices. Thus multiple edges and loops are allowed 
in graphs and digraphs. A subgraph H in a (directed) graph G is a (directed) graph 
on the same set of vertices whose set of edges is a subset of edges of G. We will 
write H C G to denote that H is a subgraph of G. For a subgraph H C G, let 
G\H denote the complement subgraph whose edge set is complementary to that of 
H . Also we will write e S G to show that e is an edge of the graph G. 

Let G be a digraph on the set of vertices 0,1, ... ,n. The vertex will be the 
root of G. The digraph G is determined by its adjacency matrix A = {cLij)o<i,j<m 
where ay is the number of edges from the vertex i to the vertex j. We will regard 
graphs as a special case of digraphs with symmetric adjacency matrix A. 

An oriented spanning tree T of the digraph G is a subgraph T C G such that 
there exists a unique directed path in T from any vertex i to the root 0. The number 
Nq of such trees is given by the Matrix-Tree Theorem, e.g., see |Sta2l Section 5.6]: 



(1) 



N G =detL G , 
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where Lg = (Uj)i<i,j<n the truncated Laplace matrix, also known as the Kirkhoff 
matrix, given by 

{^2 a ir for i=j, 
re{0,...,n}\{i} 
-a i:j for i ^ j. 

If G is a graph, i.e., A is a symmetric matrix, then oriented spanning trees defined 
above are exactly the usual spanning trees of G, which are connected subgraphs of 
G without cycles. 

For a subset I in {1, ... ,n} and a vertex i £ I, let 

di{i) =^2aij, 

m 

i.e., di(i) is the number of edges from the vertex j to a vertex outside of the subset 
I. Let us say that a sequence b — (b\, . . . , b n ) of non-negative integers is a G -parking 
junction if, for any nonempty subset / C {1, . . . , n}, there exists i £ I such that 
bi < di(i). 

If G = K n +i is the complete graph onn+1 vertices then i4T„+i-parking functions 
are the usual parking functions of size n defined in the beginning of this section. 

Theorem 2.1. cf. |Gablj The number of G -parking functions equals the number 
Nq = det Lq of oriented spanning trees of the digraph G. 

Interestingly, G-parking functions are related to the abelian sandpile model in- 
troduced by Dhar l)linT|. In the appendix we will discuss the sandpile model and 
show that Theorem 12. II is essentially equivalent to the result of Gabrielov [Gabll 
Eq. (21)] on sandpiles. In Section |H1 we will prove Theorem 12. II without using the 
sandpile model. 

We can reformulate the definition of G-parking functions in algebraic terms as 
follows. Throughout this paper we fix a field K. Let Xq = (to/) be the monomial 
ideal in the polynomial ring K[xi, . . . , x n ] generated by the monomials 

(3) m z = n^ W , 

i£l 

where / ranges over all nonempty subsets / C {1, . . . , n}. Define the algebra Ag 
as the quotient Ag = K[a;i, . . . , x n ]/Xa- 

A non-negative integer sequence b = (b±, . . . , b n ) is a G-parking function if and 
only if the monomial X - — - X • • • 3J„ IS nonvanishing in the algebra Ac- 

For a monomial ideal X, the set of all monomials that do not belong to X is a 
basis of the quotient of the polynomial ring modulo X, called the standard monomial 
basis. Thus the monomials x b , where b ranges over G-parking functions, form the 
standard monomial basis of the algebra Aq. 

Corollary 2.2. The algebra Ac is finite- dimensional as a linear space over K. Its 
dimension is equal to the number of oriented spanning trees of the digraph G: 

dim Ag — Ng- 

For an undirected graph, G-parking functions and monomials mj also appeared 
in a recent paper by Cori, Rossin, and Salvy CRS . 
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3. Power algebras 

Let G be an undirected graph on the set of vertices 0,1, ... ,n. In this case the 
dimension of the algebra Ag is equal to the number of usual spanning trees of G. 

For a nonempty subset I in {1, . . . , n}, let Dj = ^2 ieI j^j aij = J2iei ^(*) ^ c 
the total number of edges that join some vertex in / with a vertex outside of /. For 
any nonempty subset I C {1, . . . , n}, let 

(4) P/= 5>i 

\iei 

Let Jg = (pi) be the ideal in the polynomial ring K[xi, . . . , x n ] generated by the 
polynomials pi for all nonempty subsets I. Define the algebra Bg as the quotient 
B G = K[x lt . . .,x n ]/J G . 

The algebras Ag and Bg, as well as all other algebras in this paper, are graded. 
For a graded algebra A = A (B A 1 ® A 2 (B ■ ■ ■ , the Hilbert series of A is the formal 
power series in q given by 

Wi\bA = J2l k dimy^. 

fc>0 

Our first main result is the following statement. 

Theorem 3.1. The monomials x b , where b ranges over G -parking functions, form 
a linear basis of the algebra Bg- Thus the Hilbert series of the algebras Ag and 
Bq coincide termwise: Hilb^G = Hilb£>G. In particular, both these algebras are 
finite- dimensional as linear spaces over K and 

dim As = dim£>G = Ng 
is the number of spanning trees of the graph G. 

Example 3.2. Let n — 3 and let G be the graph given by 

1 

(5) G = 

2 3 

The graph G has 8 spanning trees: 





The ideals Tg and Jg arc given by 

1 G = (x\,xl,X%,x\x2,x\x%,X2xl,Xixlx 5 ) , 

J G = {x\,xl,x\, (x! +x 2 ) 3 , (xi +x 3 ) 4 , (x 2 + x 3 ) 3 , (x x +x 2 + x 3 ) 2 ) . 

The standard monomial basis of the algebra Ag is {1, x\, x 2 , x 3 ,x\, x\x 2 , x 2 x 3 , x\}. 
The corresponding G-parking functions are the exponent vectors of the basis ele- 
ments: 

(0,0,0), (1,0,0), (0,1,0), (0,0,1), (2, 0,0), (1,1,0), (0,1,1), (0,0, 2). 

We have dim Ac = dim Bq = 8 is the number of spanning trees of G, and Hilb Ag = 
WiVoBg = 1 + Zq + Aq 2 . 
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We will refine Theorem 13.11 and interpret dimensions of graded components of 
the algebras Ac and Bg in terms of certain statistics on spanning trees. Let us 
fix a linear ordering of all edges of the graph G. For a spanning tree T of G, an 
edge e G G \ T is called externally active if there exists a cycle C in the graph G 
such that e is the minimal edge of C and (C \ {e}) C T . The external activity of a 
spanning tree is the number of externally active edges. Let N G denote the number 
of spanning trees TcGof external activity k. Even though the notion of external 
activity depends on a particular choice of ordering of edges, the numbers N G are 
known to be invariant on the choice of ordering. 

Let Aq and B G be the fc-th graded components of the algebras Ag and Bg, 
correspondingly. 

Theorem 3.3. The dimensions of the k-th graded components A G and B G are 
equal to 

dimA G = dimB G = N l G ^ n ~ k , 

the number of spanning trees of G of external activity \G\ — n — k, where \G\ denotes 
the number of edges of G. 



4. Examples: tree ideals and their generalizations 

4.1. Two algebras of dimension (n + l) n_1 . Suppose that G = K n+ i is the 
complete graph onti+1 vertices. As we have already mentioned, the i£„_|_i-parkmg 
functions are the usual parking functions of size n defined in the beginning of 
Section □ 

Let I n = (mi) and J n = (pi) be the ideals in the polynomial ring K[xi, . . . , x n ] 
generated by the monomials to/ and the polynomials pi, correspondingly, given by 

mi = (x h ■ ■■x ir ) n ~ r+1 , 

pi = (x h + ■ ■ ■ + x ir y^- r+1 \ 

where in both cases / = {i±, . . . , i r } runs over all nonempty subsets of {1, . . . , n}. 
Let A n = K[xx, . . . ,x n ]/l n and B n = K[x x , . . .,x n ]/J n . 

Corollary 4.1. The graded algebras A n and B n have the same Hilbert series. They 
are finite- dimensional, as linear spaces over K. Their dimensions are equal to 

dim An = dimBn = (n + l) n_1 . 

The images of the monomials x b , where b ranges over parking functions of size n, 
form linear bases in both algebras A n and B n . 

An inversion in a tree T on the n+ 1 vertices labelled 0, . . . ,n is a pair of vertices 
labelled i and j such that i > j and the vertex i belongs to the shortest path in T 
that joins the vertex j with the root 0. 

Corollary 4.2. The dimension dim .4^ = dim£>^ of the k-th graded components 
of the algebras A n and B n is equal to 

(A) the number of parking functions b of size n such that b\ + ■ ■ ■ + b„ = k; 

(B) the number of trees on n + 1 vertices with external activity (™) — k; 

(C) the number of trees on n + 1 vertices with Q) — k inversions. 
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It is well known that the numbers (A), (B), and (C) are equal, see |Krew| . The 
inversion polynomial is defined as the sum I n (q) = X^t 9^ of lnverslons 111 T over all 
trees T on n + 1 labelled vertices. Thus the Hilbert series of the algebras A n and 
B n are equal to 

ffilb An = Hilb B n = q®I n (q- 1 ). 

4.2. Two algebras of dimension I (I + kn) 11 ^ 1 . It is possible to extend the pre- 
vious example as follows. Fix two non-negative integers k and I. Let G = K^+i 
be the complete graph on the vertices 0,1, ... ,n with the edges (i, j), i,j ^ 0, of 
multiplicity k and the edges (0, i) of multiplicity I. The i^^-parking functions are 
the non- negative integer sequences b = (b\, . . . ,b n ) such that, for r = 1, . . . , n, 

#{i\bi<l + k(r-l)}>r. 

The definition of these functions can be also formulated as Cj < I + (i — 1) fc, where 
ci < ■ ■ • < c n is the increasing rearrangement of elements of b. Such functions 
were studied by Pitman and Stanley |PiStj and then by Yan | Yan| . These authors 
demonstrated that their number equals 1(1 + fcn)" -1 . One can show, using for 
example the Matrix- Tree Theorem Q), that the number of spanning trees in the 
graph K^li equals I (l + kn)"-- 1 . Thus Theoreml2~Tl recovers the above formula for 
the number of JC^^-parking functions. 

Let I n ,h,l = (mi) and J n ,k,i — (pi) be the ideals in the ring K[xj., . . . , x n ] 
generated by the monomials to,/ and the polynomials pi, correspondingly, given by 

m I = (x il ---x ir ) l+k ^- r \ 

PI = (X h + ■ ■ ■ + Xl J«+ k ^-r)) j 

where in both cases I = {i\, . . . , i r } runs over all nonempty subsets of {1, . . . , n}. 
Let An,k,i = . . . ,x n ]/l n ,k,i and B n ,k,l = K[xi, . . . , x n ]/ J n ,k,h 

Corollary 4.3. The graded algebras A Ut k,i and B n ,k,l have the same Hilbert series. 
They are finite-dimensional, as linear spaces over K. Their dimensions are 

dim An,k,i = dhnB n ,k,i =l(J + kn) n ^ 1 . 

The images of the monomials x b , where b ranges over K^'^-parking functions, form 
linear bases in both algebras A n ,k,i and B nt k,i- 

5. Monotone monomial ideals and their deformations 

A monotone monomial family is a collection M. — {mj \ I € £} of monomials 
in the polynomial ring K[xi, . . . ,x n ] labelled by a set E of nonempty subsets in 
{1, . . . , n} that satisfies the following three conditions: 

(MM1) For I G E, to,/ is a monomial in the variables Xi, i € I. 

(MM2) For I, J 6 E such that I C J and i £ I, we have deg^m./) > deg^,. (toj). 

(MM3) For I, J € E, lcm(m/, mj) is divisible by tuk for some K D I U J in E. 

The monotone monomial ideal T = (M) associated with a monotone monomial 
family M. is the ideal in the polynomial ring K[xi, . . . , x n ] generated by the mono- 
mials mi in M.. 

It follows from (MM1) and (MM2) that condition (MM3) can be replaced by the 
condition: For /, J S E there is K D I U J in E such that tuk is a monomial in the 
Xi, i G I U J. This condition is always satisfied if I, J G E implies that I U J G E. 
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The monomial ideal Xq constructed in Section|21for a digraph G is monotone. In 
this case E is the set of all nonempty subsets in {1, . . . , n} and to/ is given by 

Remark that two different monotone monomial families may produce the same 
monotone monomial ideal. For example, the ideal Xq, for the graph G shown on J5|, 
has generator m^^y = x\x\. This generator is redundant because it is divisible 
by 2 3}= x\x^- Thus the same ideal corresponds to the monotone monomial 
family with E = {{1}, {2}, {3}, {1, 2}, {2, 3}, {1, 2, 3}}. 

Let / = {ii, . . . ,i r }. For a monomial to £ K.[xi 17 . . . , Xi r ], an I-deformation 
of to is a homogeneous polynomial p £ , . . . , Xi r ] of degree deg(p) = deg(m) 

satisfying the generosity condition 

(6) K[x il ,...,x ir ] = (R m )@(p), 

where {R m ) is the linear span of the set R m of monomials in K.[xi 1 , . . . , Xi r ] which are 
not divisible by to, (p) is the ideal in K[xj x , . . . , Xi r ] generated by p, and "©" stands 
for a direct sum of subspaces. Notice that the generosity condition is satisfied for a 
Zarisky open set of polynomials in K[xi l; . . . ,Xi T ] of degree deg(TO-). For example, 
the polynomial p = a x\ + b X2 is a {1, 2}-deformation of the monomial m = x\ if 
and only if a ^ 0. 

The following lemma describes a class of /-deformations of monomials. 

Lemma 5.1. Let I = {ii, . . . ,i r }, let m be a monomial in K[xj i; . . . , Xi r ], and let 
oti, . . . , a r £ K \ {0}. Then the polynomial 

p={a 1 x ll +--- + a r x lr f^ m 

is an I-deformation of the monomial to. 

Proof. Let m — x" 1 • • • x^ r . The generosity condition iJB} is equivalent to the con- 
dition that the operator 

A : K[x it ,...,x ir ] — > K[x it , . . . , x ir ] 
A : / h- (d/dx n )^ ■ ■ ■ (d/dx ir y>-( v . /) 

has zero kernel. Let us change the coordinates to y\ = Xi 1 ,..., y r -\ = Xi r _ 1 , 
y r = a\Xi 1 + • • ■ + a r Xi r . The operator A can be written in these coordinates as 

A(f) = (9i + ai ~d r ) ai ■ ■ ■ (flU + Or-idr)**- 1 (a r d r r(y a r 1+ -+^ ■ /) 

where dj = d/dyj. Then A(f) = c • / + g, where c is a nonzero constant and 
deg yr (g) < deg yr (f). Thus, in an appropriate basis, the operator A is given by a 
triangular matrix with nonzero diagonal elements. This implies that Ker A = 0. □ 

A deformation of a monotone monomial ideal X — {mj \ I £ E) is an ideal J = 
(pi | / £ E) generated by polynomials pj such that pj is an /-deformation of To/ for 
each / £ E. For example, according to Lemma T5. II the ideal Jq given in Sectional 
is a deformation of the monotone monomial ideal Xq- 

Theorem 5.2. LetX be a monotone monomial ideal, and R be the standard mono- 
mial basis of the algebra A = K[xi, . . . ,x n ]/X, i.e., R is the set of monomials that 
do not belong to X. Let J be a deformation of the idealX, and B = K[xi, . . . , x n ]/J. 
Then the monomials in R linearly span the algebra B. 

Remark that the set of monomials R may or may not be a basis for B. 
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Corollary 5.3. Let X be a monotone monomial ideal, J be a deformation of the 
ideal I. A = K[xi, . . . ,x n ]/1, and B — "K[x\, . . . ,x n ]/S- Then we have the follow- 
ing termwise inequalities for the Hilbert series: 

Hilb 1 < Hilb J or equivalently, Hilb A > Hilb B. 

In some cases the Hilbert series are actually equal to each other. According to 
Theorem l3.ll Hilb„4c = Hilb£>G, for any graph G. However in general the Hilbert 
series may not be equal to each other. It would be interesting to describe a general 
class of monotone monomial ideals and their deformations with equal Hilbert series. 

There is an obvious correspondence between the generators mi of a monotone 
monomial ideal T and the generators pi of its deformation J . Notice however 
that (except for very special cases) the monomial generator mj does not belong 
to the boundary of the Newton polytopes of it polynomial deformation pi. Thus 
the monomial mi is not the leading term of the polynomial pi for any term order. 
This shows that the above results cannot be tackled by the standard Grobner bases 
technique. 

6. Syzygies of order monomial ideals 

In this section we introduce a class of ideals that extends monotone monomial 
ideals and construct free resolutions for these ideals. 

Let P be a finite partially ordered set, or poset. Let A4 — {m u | u G P} be a 
collection of monomials in the polynomial ring K[a;i, . . . , x n ] labelled by elements 
of the poset P. Also let M u denote the set of all monomials divisible by m u . Let us 
say that AA is an order monomial family, and the ideal I = (A4) generated by the 
monomials m u is an order monomial ideal, if the following condition is satisfied: 

(OM) For any pair u,v G P, there exists an upper bound w G P of u and v such 
that M u n M v C M w , i.e., m w divides lcm(m tI , m v ). 

Here an upper bound means an element w such that w > u and w > v in P. In 
particular, this condition implies that the poset P has a unique maximal element. 

Every monotone monomial family is an order monomial family labelled by the 
set P = £ of subsets in {1, . . . , n} partially ordered by inclusion. Indeed, condition 
(MM3) is equivalent to condition (OM). 

Let S = K.[x±, . . . , x n ] be the polynomial ring. For a non-negative n- vector 
a = (ax, . . . , a„) £ Z™, let S(— a) = x a S denote the free Z n -graded S'-submodule in 
S generated by the monomial X — X ^ * ■ ■ X „ " . This submodule is isomorphic to S 
with the Z ra -grading shifted by the vector a. If a > b componentwise then S(—a) 
is a submodule of S(—b) and we will write S(— a) <—> S(—b) to denote the natural 
multidegree-preserving embedding of S'-modules. 

For an order monomial family A4 = {m u \ u G P} and a subset U C P of 
elements of the poset P, let mjj — lcm(m„ u G U) be the least common multiple 
of the monomials m u , u G U. We assume that m$ — 1. Also let a\j G Z n be the 
exponent vector of the monomial mu . 

Let us define the homological order complex C*(A4) for an order monomial ideal 
I = (M) as the sequence of Z n -graded S'-modules 

(7) ■•■ ^C 3 ^C 2 ^Ci^C = 
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whose fc-th component is 

C fe = (J) S(-d{ ult ... iUh y), 

I"l<"'<Mfc 

where the direct sum is over strictly increasing fc-chains u\ < ■ ■ ■ < Uk in the poset 
P. The differential dk ■ Ck — > G'k~i is defined on the component S(— as 
the alternating sum dk = 2i=i( — 1)* -^i °f the multidegree-preserving embeddings 
Ei : S(-a,{ Uu ... tUk y) S'(-a{ Ul) ... i {r i ,...,« fc }) of S-modules, where tti, . . . ,Ui, . . . , Uk 
denotes the sequence with skipped i-th element. 

Theorem 6.1. Let M. be an order monomial family. The homological order com- 
plex C*(M) is a free resolution of the order monomial ideal X = {M) . 

If m {ui,...,u k } 7^ m {ui....,Si....,u k }j for any increasing chain u\ < • • • < Uk in the 
poset P and i = 1, . . . , k, then the homological order complex C*(M) is a minimal 
free resolution of the order monomial ideal X — {M.) . 

The above construction of the homological order complex C*{M) is an instance 
of the general construction of cellular complexes for monomial ideals due to Bayer 
and Sturmfels BaSt . Their cellular complexes are associated with cell complexes, 1 
whose faces arc marked by certain monomials. In our case the cell complex is the 
geometrical order complex A = A(P) of the poset P. It is the simplical complex 
whose faces correspond to nonempty strictly increasing chains in P: 

A(P) = {{m, ■•• ,u k }CP\ui<---<u k , k > 1}. 

For example, if P consists of all nonempty subsets in {1, . . . , n} ordered by inclusion 
then A(P) is the barycentric subdivision of the (n — l)-dimensional simplex. The 
face F of A(P) given by an increasing chain u± < ■ ■ ■ < Uk is marked by the 
monomial mp = m/ Ul UJt \. For a monomial m, let A< m denote the subcomplex 
of A(P) formed by the faces F whose mark mp divides m: 

A< m = {F e A(P) | m F divides to}. 

The faces of A(P) are partially ordered by containment of closures. More precisely, 
F > F' if the increasing chain for the face F' is a subchain of the chain for F. 

A result of Bayer and Sturmfels [BaStl Proposition 1.2] on cellular complexes 
implies the following statement. 

Lemma 6.2. The complex C*(A4) is exact if and only if A< m is acyclic over K 
for any monomial to. //, in addition, ttif ^ mpi , for any F > F' , then the complex 
C<t{M) is a minimal free resolution of the ideal X = (M) ■ 

Actually, the subcomplex A< m is not only acyclic but also contractible. This 
follows from the following result of Narushima. 

Lemma 6.3. Nary] Let M u , u £ P, be a finite collection of subsets in some set 
971, whose index set is a poset P. Assume that, for any u,v £ P, M u n M v C M w 
for some upper bound w £ P of u and v. Then, for any to £ 971, the subcomplex 
of the order complex A(P) of P formed by the following collection of nonempty 
increasing chains in P 

{{ Ul < ■ ■ ■ < u k } | m g M U1 n • ■ • n M Uk , k > 1} 

1 "Cellular complexes" should not be confused with "cell complexes." The former are homo- 
logical complexes and the latter are geometrical complexes. 
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is contractible. 

Theorem 16. II easily follow from Lemmas 16 . 21 and 16 . 31 

Proof of Theorem \b\l\ Let M = {m u \ u £ P} be an order monomial family. In 
view of Lemma l6.2l it is enough to show that the subcomplex A< m is contractible for 
any monomial m. According to (OM), the conditions of Lemma f6. 31 are satisfied, 
where 9Jt is the set of all monomials in K[xi, . . . , x n \. For an increasing chain 
U\ < • • • < lit in P, the intersection M Ux PI • • • H M Uk is the set of monomials 
divisible by rn^ Ult Uk y. Thus the contractible complex from Lemma 16.31 is exactly 
the complex A< m . Lemma l6~2l implies that the homological order complex C*(A4) 
is exact. □ 

Let us now assume that M. — {mj \ I G £}, mj — YiiPi x i T ^ > ^ s a monotone 
monomial family. In this case, for a strictly increasing chain of subsets 1\ C • ■ ■ C I k , 
the least common multiple msj^ j k \ = lcm(m/ 1 , . . . , mj k ) is given by 

(8) "Kw^n^" n < j2(i2) x---x n <y ik) - 

h£h I2£l2\li i k Gl k \I k -i 

In other words, the exponent vector CLii lt ,,, t i h \ = ( a ij • • • j a n) € Z n of the monomial 
™-{Ji,.. is given by 

( vi r (i) if i G I r \ I r -i, 
a% = < 

I if »07*, 

where we assume that Iq = 0. 

Let us say that the monotone monomial family M. and the corresponding mono- 
tone monomial ideal X = (M.) are strictly monotone if the following additional 
conditions hold: 

(SMI) The ideal X is minimally generated by the set of monomials {mi \ I G £}, 

i.e., there are no two elements 7 ^ J in £ such that m/ divides mj. 
(SM2) For any 1 C J C if in E, there exists i E J \ I such that > vk{^)- 

For example, a monotone monomial family such that the inequality in (MM2) is 
always strict and fi{i) > 0, for any / G E and i G 7, will be strictly monotone. 

Conditions (SMI) and (SM2) are equivalent to the statement that m/^ j^i 7^ 
TO {/i 1 i k } f° r anv increasing chain 7i C • • • C 1^ in E and £ = 1, . . . , ft. Thus 
Theorem 16. II specializes to the following statement. 

Corollary 6.4. Lei M te a monotone monomial family. Then the homological 
order complex C*(A4) is a free resolution of the ideal X — {M). If M is a strictly 
monotone monomial family, then C*(M) is a minimal free resolution of the strictly 
monotone ideal X. 

Homological order complexes are related to Scarf complexes of generic mono- 
mial ideals. Let X = (mi, . . . , m r ) be an arbitrary ideal in the polynomial ring 
S = K[xi, . . . , x n ] minimally generated by monomials mi, . . . , m r , and let mjj = 
lcm(m u I u G U) for U C {l,...,r}. The geometrical Scarf complex of X was 
defined by Bayer, Peeva, and Sturmfels BPSj as the following simplicial complex: 

A scari = {u C{l,...,r}\m v ^m v for all subsets V ^ U, and \U\ > 1}. 

The corresponding cellular complex is called the homological Scarf complex. 
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Lemma 6.5. Let X = (m u \ u G P) be an order monomial ideal. Then the geomet- 
rical Scarf complex Aj Car ^ is a subcomplex of the geometrical order complex A(P). 

Proof. Let U be a subset of elements in P. Suppose that U contains two incom- 
parable elements u and v. Let us pick an upper bound w of u and v provided by 
condition (OM). Let V = U U {w} if w £ U, or V = U \ {w} if w G U. Then, 
according to (OM), mjj = my. Thus V does not belong to the geometrical Scarf 
complex. 

This implies that, for any U G Aj , all elements of U are comparable with 
each other, i.e., U is an increasing chain in the poset P. Thus U G A(P). □ 

Let us say that a monomial m strictly divides a monomial to', if to divides to' 
and deg x . (ml /m) ^ whenever deg x .(m') 7^ 0. According to Miller, Sturmfcls, 
and Yanagawa [MSY| . the ideal X = (mi, . . . , m r ) is called generic monomial ideal 
if the following condition holds: 

(GM) If two distinct minimal generators m u and m v have the same positive degree 
in some variable Xi, there is a third generator m w which strictly divides 
lcm(TO u , m v ). 

The general, the Scarf complex may not be acyclic, but, for generic monomial 
ideals, Miller, Sturmfels, and Yanagawa |MSY| proved the following result. 

Proposition 6.6. jMSYI Corollary 1.8] If I is a generic monomial ideal, then the 
homological Scarf complex is a minimal free resolution ofX. 

We will see in Section that there are strictly monotone monomial ideals that 
are not generic and there are generic monomial ideals that are not strictly mono- 
tone. The following claim shows that these two classes of ideals have an interesting 
intersection. 

Proposition 6.7. Let X — (mi \ I G E) be a monotone monomial ideal such that 
the inequality in (MM2) is always strict, and vi(i) > 0, for any I G S and i G I. 
Then the monomial ideal X is both generic and strictly monotone. In this case, the 
geometrical/homological order complex for X coincides with geometrical/homological 
Scarf complex for X. 

Proof. If monomials m/ and toj have the same positive degree in some variable 
then I and J are incomparable in E: I % J and J % I. By (MM3) there exists 
K D I U J such that m^ divides lcm(TO/, toj). Then K ^ I, J. Since we assume 
that the inequality in (MM2) is strict, mx strictly divides \cm(mi,mj). It follows 
that the ideal X is generic. 

According to Lemma Ifi. 51 the geometrical Scarf complex Aj arf is a subcomplex 
of the geometrical order complex A(E). Let us prove that Aj carf = A(E). We need 
to show that, for any increasing chain I\ C • • • C in E, we have m^j lt j k j ^ mu, 
where R is any subset of E such that R ^ {I\ , . . . , Ik}- This is clear if R is a subchain 
in ii C • ■ ■ C If.. Otherwise, suppose that m{j lr .. i j Ji j = m^ and R contains an 
element J £ {ii,...,/^}. Then toj divides l rn{i 1 ,...j k }- According to conditions 
of the proposition, the monomial mj depends nontrivially on all x%, i G J. Thus 
J C If.. Then J C I r and J $Z I r -i for some r G {1, . . . , k}. (We assume that 
Iq = 0.) Pick an element i G J\I r -i- Then vj(i) < deg x .(m^ Il ^ j k j) = vj r (i) 
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because mj divides m^] lt j k y. Also we have vj(i) > vi r (i) because J C I r . 
Contradiction. □ 

The k-th Betti number f3 k {X) of an ideal X is the rank of the fc-th term in a 
minimal free resolution of X. The graded Betti number f3 k ,d(2) of a graded ideal X 
is the number of direct summands in the fc-th term of a minimal free resolution of 
X with generator of degree d. Then (3 k (X) = J^d Pk,d@)- 

Let d{I\ 1 . . . , I k ) be the degree of the monomial j fc \ given by 

(9) d(i u ...,i k )=j2 v h(h)+ ^ 2 (*2) + ---+ vi ^- 

hell i2&h\h i k el k \lk-i 

Corollary 6.8. Let X = (mi I £ E) be a strictly monotone monomial ideal. The 
k-th Betti number (3k (X) ofX is equal to the number of strictly increasing k-chains in 
the poset E. Moreover, the graded Betti number /3 k> d(I) °fZ * s equal to the number 
of strictly increasing k-chains I\ C • ■ ■ C I k in E such that d(I\, . . . . Ik) = d. 
In particular, if E is the set of all nonempty subsets in {1, . . . , n} then 

(10) /3 k (I) = MS n+1 , k+1 , 

where SVi+i.fc+i is the Stirling number of the second kind, i.e., the number of parti- 
tions of the set {0, . . . , n} into k + 1 nonempty blocks. 

The last claim is obtained by counting strictly increasing fc-chains of nonempty 
subsets I\ C • • • C Ik in {1, . . . ,n}. Indeed, these chains are in one-to-one cor- 
respondence with partitions of the set {0, . . . ,n} into a linearly ordered family of 
k + 1 nonempty blocks (ii, I2 \ I\, • ■ • ,Ik \ Ik-i, {0, . . . ,n} \ Ik) such that the last 
block contains 0. There are fc! ways to pick such a linear ordering of blocks. 

Let us say that a (directed) graph is saturated if all off-diagonal entries of the 
adjacency matrix A = (a^) are nonzero: dy > 1 for i ^ j. If G is a saturated 
digraph then the monotone monomial ideal Xq constructed in Section [5] satisfies 
the conditions of Proposition 16. 71 

Corollary 6.9. The monotone monomial ideal Xq, for a saturated digraph G, is 
both strictly monotone and generic. In this case E is the poset of all nonempty 
subsets in {1, ... ,71}. The homological order complex C*(M.), which coincides with 
the homological Scarf complex, gives a minimal free resolution of the ideal Xq. Its 
Betti numbers are given by formula <|10t • 

It would be interesting to find a minimal free resolution of the ideal Xq for any 
non-saturated digraph G. More generally, it would be interesting to find a minimal 
free resolution for any monotone monomial ideal. 

Computer experiments suggest the following conjecture on Betti numbers of 
deformations of monotone monomial ideals. 

Conjecture 6.10. Let J be a deformation of a monotone monomial ideal X such 
that dimjcKfa;!, . . . ,x n ]/X — dimK EC[:ri, . . . ,x n ]/J. Then all graded Betti numbers 
of the ideals X and J coincide: (3k,d{I) — Pk.d{<J)- In particular, for a graph G, 
the ideals Xq and Jq have the same graded Betti numbers. 
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7. Examples 

Let us illustrate Corollaries 16.41 16.81 and 16.91 and Proposition 16.71 by several 
examples. In all examples n — 3, S = K.[xi, X2, X3], and S(—d) denotes the Z- 
graded S'-module isomorphic to S with grading shifted by integer d, so that the 
generator has degree d. 

Example 7.1. Let G = K4 be the complete graph on 4 vertices. This graph is 
saturated. Thus the monomial ideal Xq is both strictly monotone and generic. 
The poset S consists of all nonempty subsets in {1, 2, 3}. The Hasse diagram of S 
marked by the monomials mj is given by 



X1X2X3 




The poset £ has seven 1-chains, twelve 2-chains, and six 3-chains. In this case the 
geometrical order complex A = A(£) is the barycentric subdivision of a triangle. 
The following figure shows the complex A with faces marked by vectors a/jj i k \' 



A 



230 032 

/ 131 \ 
220 /231 I 132\022 

/ 221 111 122 \ 
320 321^^123 ^ 
311 212 113 
1 



302 



■ 203 



The Bctti numbers (/?o, Pi, 02, 0z) = (L 7, 12, 6) of the ideal X = Xk 4 , which are also 
the /-numbers of the order complex A, can be expressed in terms of the Stirling 
numbers by formula (|10fl . The graded Betti numbers of this ideal are indicated on 
the following minimal free resolution: 







S(-6) 6 — ► 5(-5) : 



5(-3) 4 ®5(-4) s 



S 



S/I 



0. 



This resolution is the homological order complex and also the homological Scarf 
complex of X. 

Similarly, a minimal free resolution of the ideal X n = Xx n+1 associated with 
the complete graph K n+ i is given by the cellular complex corresponding to the 
simplicial complex A(E) = A^ , which is the barycentric subdivision of the (n— 1)- 
dimensional simplex, cf. |MSYI Example 1.2]. 

Example 7.2. Let G be the graph given by 

1 



G 



2 3 

This graph is not saturated and the monotone monomial family that generates the 
ideal I = Xq = (xf , x%, x%, x\x^, x\x\, x^x\, x\x%x-$) will not be strictly monotone 
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if we assume that the labelling set £ consists of all nonempty subsets in {1,2,3}. 
As we mentioned before, the generator m{ 13 } = x\x\ is redundant. The same 
monomial ideal X is minimally generated by the strictly monotone monomial family 
{mi | J e £} with S = {{1}, {2}, {3}, {1, 2}, {2, 3}, {1, 2, 3}}. The Hasse diagram 
of this poset £ marked by the monomials rax is given by 



2:1X2X3 




2:3 



This poset has six 1-chains, nine 2-chains, and four 3-chains. Its geometrical order 
complex A = A(£) with faces marked by vectors ax u .. x h is shown on the following 
figure: 



(120 




300 003 



The homological order complex give a minimal free resolution of the ideal X: 

— ► S(-5) 4 — ► S'(-4) 9 — ► S(-2) 2 S(-3) 4 — ► S — >S/X — ► 0. 

In this case the ideal X is also generic and the above resolution is the homological 
Scarf complex. 

Example 7.3. Let X — (a^x%,x\xs,x\X<2,x%) be the monotone monomial ideal, 
whose poset £ marked by the monomials is given by 



X1X2X3 



S = 




2 2 2 

The ideal X is strictly monotone but is not generic. The geometrical order complex 
with faces marked by vectors asx x ,...,i k } is given by 

, 220 221 111 121 021 
A = • • • 

It produces the following minimal free resolution of the ideal X: 

— ► 5(-4) © 5(-5) — > S(-3) 2 © S(-4) — ► S — ► S/X — ► 0. 

On the other hand, the geometrical Scarf complex in this case is disconnected and 
docs not give a resolution for X. 
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Example 7.4. Let G be the graph given by 

1 



G = 




2 3 

Then X = Xq — (x\ , x\, x\, x\x\, x\x-z, x\ x\, X1X2X3Y The Hasse diagram of S 
marked by the monomials rrij is given by 

X1X2X3 



s = 




^,3 ~,2 

J, 1 J-2 Xg 



The corresponding homological order complex C*(Ai) gives a free resolution, which 
is not minimal. In this case the monomial generator msin = x\x\ is redundant 
because it is divisible by mti} — x\. However, the family M \ {mu^x} is not a 
monotone monomial family because condition (MM3) no longer holds. On the other 
hand, X is a generic ideal and its Scarf complex gives a minimal free resolution: 

— ► S"(-5) 4 — ► S(-4) 9 — ► S(-2) 2 ® S(-3) 4 — ► S — ► S/X — ► 0. 



In Examples 17 . II [7"2l and 17.41 above, the graded Betti numbers of the deformed 
ideal Jq coincide with graded Betti numbers of Xq. 



8. HlLBERT SERIES AND DIMENSIONS OF MONOTONE MONOMIAL IDEALS 

In this section we give formulas for the Hilbert series and dimensions of monotone 
monomial ideals. Then we prove Theorem 12. II 

Let {mi | / G E}, to/ = Yii^i^'^^ be a monotone monomial family, let 
2 = (mi) be the corresponding ideal in the polynomial ring K[xi, . . . , x n ], and let 
A = K[xi, . . . ,x n ]/I. 

Proposition 8.1. The Hilbert series of the algebra A equals 





Hilb.4 



i + ^(-i) fc E q d(Il - Ik) 

k>l liQ-CI k 



(l-g)« 

where the sum is over all strictly increasing chains in E and the number d(I\ , . . . , Ik) 
is given by 10- 

First Proof of Provosition \H . 1\ According to Corollary 16.41 the homological order 
complex {7J is a free resolution of X. Thus 

Hilb^ = ^(-l) fe HilbC fe = ]T(-l) fc ]T HilbS(-a {Jl ,..., Jfc} ). 
Since Hilb S(-a {Iu ...j k} ) = j (\ _ the proposition follows. □ 
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Let us give another more expanded proof of Proposition 18.11 We will need the 
improved inclusion-exclusion formula due to Narushima |Naru| . For a subset M in 
some set 971, let x(M) denote the characteristic function of M: 



X (M) 



if a G M, 

if a G £01 \ M. 



Lemma 8.2. Narij] Let M u , u G P , be a finite collection of subsets in some set 
971, whose index set is a poset P, such that, for any u,v £ P, M u n M v C M w for 
some upper bound w of u and v. Then we have 

xW\ |J M„)= x (£SJt) + ^(-i)* An-niJ, 

u£P fe>l ui<-<n t 

where the second sum is over all strictly increasing chains u\ < ■ ■ • < itfe in the 
poset P. 

Proof. According to the usual inclusion- exclusion principle, see {Stall Section 2.1], 
we have 

x(m \ |J m u ) = x (m) - £ x(Mu) + n M v ) — . 

kGP u u,v 

The general summand in this expression is sjj — (— l) k x(M Ul fl • • • fl M Uk ), where 
U = {u\, . . . , Ufc} is an unordered fc-element subset in P. We argue that if we take 
the summation only over increasing chains u\ < ■ ■ ■ < Uk in P we get exactly the 
same answer. Indeed, let us show that the contribution of all other subsets U is 
zero. We will use the involution principle, see {Stall Section 2.6]. Let us construct 
an involution i on the set of all subsets U C P of all possible sizes k > such that 
the elements of U cannot be arranged in an increasing chain. Let us fix a linear 
order on elements of P. Find the lexicographically minimal pair of of incomparable 
elements u and v in U, i.e., u ^ v and v ^ u. Let id £ F be the minimal (with 
respect to the linear order) upper bound of u and v such that M u fl M v C M w . 
Define the map t as follows: 

f UU{w}, if w£U, 
\ U\{w}, if w G U. 

Then i is an involution such that |t(i7)| — \U\ ± 1. Conditions of the lemma imply 
that sjj = — 3 t (m. Thus all summands su corresponding to non-chains cancel each 
other. □ 



U i 



Second Proof of Provosition \H.l\ Let 9Jt be the set of monomials in K[xi, . . . , x> 
and, for I G S, let Mj C 971 denote the set of monomials in K[xi, . . . , x n ] divisible 
by raj. For a subset of monomials M C Wl, let 

[M] = J2 1 dCe(m) = J2 q dce{m) x(M){m). 

m£M mSOT 

Then Hilb^l = [M \ \JM{\ and [971] = (1 - q)- n . All conditions of Lemma IO 
are satisfied, where P = S. For an increasing chain 1\ C ••• C 1^, the least 
common multiple (JSj of the monomials , . . . , mi k has degree d{I\ , . . . , /&) . Thus 
[M h n- ■ -HAf/J = PropositionlOfollows from LemmalO □ 
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Lemma 8.3. The algebra A is finite- dimensional as a linear space over K if and 
only if E contains all one-element subsets in {1, ... ,n}. 

Proof. If there is i € {1, . . . , n} such that {i} $ £ then the powers x\ form an infinite 
linearly independent subset in A. Thus A is infinite-dimensional. Otherwise, if £ 
contains all one-element subsets, then the algebra A is finite-dimensional. Indeed, 
a monomial x^ 1 . . . x^™ vanishes in A unless a\ < kmi (1), . . . ,a n < vt n \(n). □ 

Proposition 8.4. j4sswme £/ia£ £ contains all one-element subsets in {1, . . . ,n} 7 
and let v(i) = vu\[i). The dimension of the algebra A is given by the following 
polynomial in the variables | / £ S,i G /}: 

dirndl = Y, (-!)* II (Kii)-^i(ii)) x II 

(11) 

x---x J| (y(i fc ) -v Ik {i k )) x f(ife+i), 
i k eik\ik-i ik+i&h 

where the sum is over all strictly increasing chains I\ ^ • • • C of nonempty 
subsets in {1, . . . , n} of all sizes k > 0, including the empty chain of size k = 0. 

Proof. Let 9Jt be the set of monomials x" 1 ■ ■ ■ x^™ such that a\ < ^(1), ... ,a n < 
v(n). A monomial x a vanishes in the algebra A unless x a € OJl. Let M/ = MjD^Xft. 
Lemma f8. 21 for 971 implies that 

dim^= | art | + ^(-i) fe Y l^n-nM/j. 

k>l /iC-.-C/ fc 

The intersection D • • • H M/ fc is the set of all monomials in VJl divisible by 
the monomial mi lt ...j k given by 1(8}. Thus (— l) fe \Mj x f~l • • • PI Mj k \ is equal to the 
summand in 1|11|) . □ 

Remark that if 1%, . . . , is not a chain of subsets then \Mj 1 f) ■ ■ ■ n Mj fc | may 
not be a polynomial in the vi(i). It may include expressions like min(z/j(i), 
Thus the inclusion-exclusion principle does not immediately produce a polynomial 
expression for dim .4. Miraculously, all non-polynomial terms cancel each other. 

We can now prove Theorem 12.11 that claims that dimension of the algebra Ag 
equals the number of oriented spanning trees of G. 

Proof of Theorem \2.1\ cf. Gabl , Appendix E] Let G be a digraph on the vertices 
0, . . . ,n, and let A = (a y -) be its adjacency matrix. Specializing Proposition 18.41 
we obtain the following polynomial formula for the dimension of the algebra Ag '■ 



(12) 



dim^l G = J2 II J2 "iih x II E a ^h 

x • • • x n \J2 a ^ x n n a ik+ljk+1 ] , 

«*e/ h \/fc_i \3fc6Ijb / i fe+1 e{l,...,n}\4 \0<] k+ i<n 



where the sum is over all strictly increasing chains I± C • ■ ■ C 4 of nonempty 
subsets in {1, ... ,n} of all sizes k > 0. In this formula, we assume that an = 0. 
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Let us show that the expression (|12(l for dim^lG is equal to the number of 
oriented spanning trees of G. We will use the involution principle again. 

Let us first give a combinatorial interpretation of the right-hand side of (|12fl . 
The summand that corresponds to an increasing chain 1* = I± C • • • C I/. is equal 
to (— 1) times the number of subgraphs H of G such that 

(a) H contains exactly n directed edges (i, /(«)) for i = 1, . . . , n. 

(b) If i £ I r \ I r -i then f(i) £ I r . (We assume that I = 0.) 

For such a subgraph H, let Jh Q {1, . . .,n} be the set of vertices i such that 
f p (i) = for some power p, i.e., Jh is the set of vertices i such that there is a 
directed path in H from i to the root 0. Notice that if i £ {Jl r then f p (i) £ IJir-, 
thus f p (i) ^ 0, for any p. Thus Ji, . . . ,1/. C J# = {1, . . . , n} \ J H . Also notice that 
H is an oriented spanning tree of G if and only if Jh = {1, . . . , n}. 

Let us now construct an involution n on the set of pairs (/* , H ) such that -ff is 
not an oriented spanning tree. In other words, the involution n acts on the set of 
pairs (I* , H ) with nonempty Jh ■ It is given by 

f (/iC-CJ H ,fl) if/fc=7ff, 

1 {hC...CI k CJ H ,H) if/fc^Jff. 

The contribution of the pair (/* , H ) to the right-hand side of l|12fl is opposite to the 
contribution of H). Thus the contributions of all subgraphs H which are not 
oriented spanning trees cancel each other. This implies that dim .4(3 is the number 
of oriented spanning trees. □ 

It would be interesting to find a combinatorial proof of Theorem 12. II In other 
words, one would like to present a bijection between G-parking functions and 
oriented spanning trees of G. There are several known bijections between the 
usual parking functions and trees. One such bijection is relatively easy to con- 
struct. There is a more elaborate bijection that maps parking functions b with 
bi H h b n = k to trees with (™) — k inversions, see |Krewj . 



9. Square-free algebra 

Let G be a graph on the set of vertices 0, . . . , n. We will say that a subgraph 
H C G of the graph G is slim if the complement subgraph G \ H is connected. Let 
us associate commutative variables <j> e , e £ G, with edges of the graph G, and let 
$G be the algebra over IK generated by the <fi e with the defining relations: 

(4> e ) 2 = 0, for any edge e; 

Y[ 4>e — 0, for any non-slim subgraph H C G. 

eeH 

Clearly, the square- free monomials 4>h — Ileeff wnere H ranges over all slim 
subgraphs in G, form a linear basis of the algebra $q. Thus the dimension of $g 
is equal to the number of connected subgraphs in G. 
For i = 1, . . . , n, let 

Xi = C ly£ 4> e , 

eGG 
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where 

{1 if e = (i,j), i < j; 
-1 if e = (i,j), i > j; 
otherwise. 

Define Cq as the subalgebra in <&q generated by the elements X\ , . . . , X n . 

Fix a linear ordering of edges of the graph G. Recall that N G denotes the number 
of spanning trees of G with external activity k, see Section [21 

Theorem 9.1. (1) The dimension of the algebra Cq as a linear space over IK 
equals the number of spanning trees in the graph G. 

(2) The dimension of the k-th graded component C G of the algebra Cq equals the 
number Nq ™ k of spanning trees of G with external activity \G\ — n — k. 

Recall that, for a nonempty subset / C {1, . . . , n}, Dj = J2iei j& ai i 1S ^ e 
number of edges in G that connect a vertex inside I with a vertex outside of /, see 
Section 

Lemma 9.2. For any nonempty subset I C {1, . . . , n\, the following relation holds 
in the algebra Cq: 




This lemma shows that the algebra Cq is a quotient of the algebra Be- We will 
eventually see that Bq — Cq, but let us pretend that we do not know this yet. 

Proof. Let Hi C G be the subgraph that consists of all edges that connect a 
vertex in / with a vertex outside of /. We have Y^iei -^-i = J2 e eH T ^fie- Thus 
(Xhsi Xi) D ' = ± rieG-ff/ ( t >e ~ 0, because Hi is not a slim subgraph of G. □ 
Let Sg be the subspace in K[j/i, . . . , y n ] spanned by the elements 

dH = \\ a e , 

for all slim subgraphs H C G, where a e — yi — yj, for an edge e = (i, j), i < j. Let 
S G denote the fc-th graded component of the space Sg- 

Lemma 9.3. For any graph G and any k, we have dimC G = dimiS^. 

Proof. Let bn, a be the coefficient of J| e6 H 4> e in the expansion of X" 1 ■ • ■ , where 
a = (ax, . . . , a„). Then dimC G is equal to the rank of the matrix B = (bn,a), where 
H ranges over all slim subgraphs in G with k edges and a ranges over non-negative 
integer n-element sequences with ai + ■ ■ ■ + a n = k. On the other hand, ba, a is also 
equal to the coefficient of y" 1 • • • in the expansion of an- Thus dimS G equals 
the rank of the same matrix B = (6jj )0 ). □ 

For a spanning tree T in G, let T + denote the graph obtained from T by adding 
all externally active edges. In virtue of Lemma 19.31 the following claim implies 
Theorem Q 

Proposition 9.4. The collection of elements &q\j-+ > where T ranges over all span- 
ning trees of G, forms a linear basis of the space Sq ■ 
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Let us first prove a weaker version of this claim. 

Lemma 9.5. The elements «g\t+> where T ranges over all spanning trees of G, 
spans the space S G . 

Proof. Suppose not. Let H be the lexicographically maximal slim subgraph of G 
such that an cannot be expressed as a linear combination of the ce G \r+ ■ Then 
there exists a cycle C = {ei, . . . , e{\ C G with the minimal clement e\ such that 
H n C = {ei}. Then a ei is a linear combination of a e2 , . . . ,a ei . Let Hi be the 
graph obtained from H by replacing the edge ei with a. For i = 2, . . . , I, the graph 
Hi is a slim subgraph of G, which is lexicographically greater than Hi . Then ctu 
can be expressed as a linear combination of an 2 , ■ • ■ > o i H l ■ Contradiction. □ 

Proof of Proposition fP^j Recall that N G denote the number of spanning trees in 
the graph G. In view of Lemma f9. 51 it is enough to show that diniiS G = N G . We 
will prove this statement by induction on the number of edges in G. 

If G is a disconnected graph then it has no slim subgraphs and dim S G = N G = 0. 
If G is a tree then dim S G = N G = 1. This establishes the base of induction. 

Suppose that G is a graph with at least one edge. Pick an edge e of G. Let G\e 
be the graph obtained from G by removing the edge e, also let G/e be the graph 
obtained from G by contracting the edge e. Then N G — N G \ e + N G / e . Indeed, 
for a spanning tree T in G, we have either e T or e G T. The former trees are 
exactly the spanning trees of G\e. The later trees are in a bijective correspondence 
with spanning trees of G/e. This correspondence is given by contracting the edge 
e. Assume by induction that the statement is true for both graphs G\e and G/e. 

Let S' G C Sg be the span of the as' with slim subgraphs H ' C G such that 
e £ H' and let <S£j C Sq be the span of the an" with slim subgraphs H" C G such 
that e ff". Then the space Sq is spanned by S' G and 5^. Thus 

(13) dim S G = dim S' G + dim - dim^ fl S G ). 

We have S' G = (yi — Hj)S G \ e , where e = Thus dim S' G = dimS G \ e . Let 

p : f(yi, ...,y„) i-> /(yi, . . . , y„) mod (y 4 - ^) be the natural projection. Then 
p(S£) = S G/e and 5^ n S' G C Ker(p). Thus 

(14) dim^c = dim5 G/e + dimKer(p) > dim5 G/e + dim(5 G n S G ). 
Combining (fT^|l and l(Ti]l. we get 

dim<S G > S G \ e + dim5 G/e . 

By the induction hypothesis, the right-hand side of this expression equals N G \ e + 
Na/e = Nq. Thus dim5 G > A^ G . On the other hand, Lemma 19.51 implies that 
dim<S G < N G . Thus dim5 G = A^ G , as needed. □ 

10. Proof of Theorems EHl EOl and 1531 

Let {mi | I £ S} be a monotone monomial family, and let X — {mi \ I S E) be 
the corresponding monomial ideal in K[xi, . . . , x n \. 

For a subset / = {ii,...,i r } S {1, . . . , n}, let K[xi] = K[xi 17 . . . , Xi k ], and let Tti 
denote the set of all monomials in the variables Xi, i E I. Also let 971 = 071(1,. ..,«.}• 
For / e S, let A// = to/ • 97? be the set of all monomials in K[xi, . . . , x n ] divisible 
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by mi. The standard monomial basis R of the algebra A = K[xi, . . . , x n ]/X is the 
set of monomials 

R = m\\J Mi 

/es 

that survive in the algebra A. 

For I, J £ S, denote by mj/j the monomial obtained from mj by removing 
all Xi's with i £ /, and let Mm = toj/j • 9Jtj, where I = {1, . . . , n} \ I. Let Xj 
be the monomial ideal in the polynomial ring K[a;y] generated by the monomials 
{ m J/i J £ J I}- It follows from the monotonicity condition (MM2) that the 
ideal 1/ is also generated by the set of monomials {mj/j | JeE,JD /}. Notice 
that Ti is also a monotone monomial ideal. Let Ri be the standard monomial basis 
of the algebra Ai = K[xj]/Xi: 

Ri = Mj\ [jMj/i. 

JDI 

Proposition 10.1. The polynomial ring K[xi, . . . , x n ] decomposes into the direct 
sum of subspaces: 

] = (R) © mi K[xi] {Ri} , 
/es 

where (R) and (Ri) denote the linear spans of monomials in R and Ri, respectively. 

Lemma 10.2. For any monomial x a = x" 1 • ■ • x ^ 1 in the ideal X there is a unique 
maximal by inclusion subset J £ £ such that x a € Mi. 

Proof. Let £' = {/ e E | x a e Mi}. If I, J € then, according to condition 
(MM3), there is an upper bound K G S of / and J such that M/ n Mj C M K . 
Thus Mk € This implies that £' has a unique maximal element. □ 

Proof of Provosition MO . 1\ For /eS, let M™ ax be the following set of monomials: 

Mf ax = Mj \ (J Mj, 

i.e., M™ 2 " is the set of monomials x a € 3JI such that / is the maximal by inclusion 
subset J 6 E with x a £ Mj, see Lemma ll().2l Thus the set of all monomials in 
K[xi, . . . , x n ] decomposes into the disjoint union 

(15) Tt = R U |j Mf ax . 

/es 

Using monotonicity condition (MM2), we obtain, for JC J, 

Mi \ Mj = mi xWIi x (9% \ M, 7// ), 

where the notation "x" means that every monomial in the left-hand side decom- 
poses uniquely into the product of monomials. Thus we have 

(16) Mf ax = P| (Mj \ Mj) =mi xTIi x f] (9% \ M JfI ) = mi x SOT/ x i? 7 . 

JDI JDI 

Formulas 1)15(1 and (|16|l imply the required statement. □ 

Let pi, I G S, be a collection of polynomials such that pi is an /-deformation of 
the monomial mi. Remarkably, a similar statement is valid for the polynomials pi. 
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Proposition 10.3. The polynomial ring K[xi, . . . , x n ] decomposes into the direct 
sum of subspaces: 

Ji\|Xl , . . . , x n } = (R)®Q)piK[xi] (Ri). 
/es 

Proposition 110.31 immediately implies Theorem 15.21 which says that the mono- 
mials in R linearly span the algebra B = K[xi, . . . , x n ]/ {pi / £ E). 

Lemma 10.4. Suppose that a polynomial p £ is an I-deformation of a 

monomial to £ K[a;j], see JBJ. Then for any polynomial f £ K[a;j] there exists 
a unique polynomial g £ K.[xi] such that the difference m ■ f — p ■ g contains no 
monomials divisible by to. The map f i— > g is one-to-one. 

Proof. According to the generosity condition JBJ the polynomial m • /, as well as 
any other polynomial in K[acj] , can be written uniquely in the form m ■ f = p- g + r, 
where g £ K[xj] and r is in the linear span (R m ) of monomials in K[x/] not divisible 
by to. This defines the map / i— ► g. 

On the other hand, for any g £ K[aij] there exist unique / G WL[xi] and r £ (R m ) 
such that p ■ g — to ■ f — r. Thus the map / i— > g is invertible. The statement of the 
lemma follows. □ 

Proof of Provosition \1(). , f fl Pick any linear ordering 1%, . . . , 7jy of the set £ compat- 
ible with the inclusion relation, i.e, the inclusion I s C It implies that s < t. Let 
£( s ) = {ii, . . . , I s } and X( s ) ~ {/ s , . . . , 1^} be initial and terminal intervals of E. 

We will prove by induction on N — s that the polynomial ring Kfaii, . . . , x n ] 
decomposes into the direct sum of subspaces 

(17) ] = (R) © mv K[x r ] (R P ) © Pi" K[xp>] (Ri») . 

i'es w /// eS < s +i) 

If s = N then l|17|l is true according to Proposition IIP. l( This gives the base of 
induction. 

Assume by the induction hypothesis that (|17|) holds for some s and derive the 
same statement for s — 1. Let I = I s . For a polynomial </> £ K[ii, . . . , x n ], write its 
unique presentation 

(18) (j) = r + ^ m i' ' fl' ' r i' + m i ■ fi - r i + Pl " ' ' ri " ' 

/'es s -i i"e^ s+1 

where r £ (R) and fj £ K[xj] and rj £ Rj, for any J £ S. 

Let // £ IK[a;/] be the unique polynomial, provided by Lemma |l 0.41 such that 
the difference d = to/ ■ fi — pi ■ f £ K[xi] contains no monomials divisible by toj. 
Let tp £ K[xx, . . . , x n ] be the polynomial obtained from </> by keeping all terms 
in (|18fl except for to/ • // • r/ which we substitute by the term pi ■ fi ■ rj . Then 

— = d ■ Tj. Pick any monomial e in d. Remind that, according to lit) II . M™ ax 
is the set of all monomials in m/K[xj] (Rj). If e • 77 ^ Af™ ax for all Jg E, then 
e ■ rj £ (R). Otherwise, suppose that e ■ 77 G M™ ax for some J. If J G! I, then 
e ■ 77 G Af™ ax C Mj implies that 77 G Mj/j, which is impossible. Thus J C I. Also 
J ^ I because e is not divisible by m/. This shows that 

- V> e (i?) © (Mr x ) . 

JC/ 
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Therefore, <j) can be written as 

(19) (p = f+ ^ m I' ' fl' 'fl'+Pl'fl-ri+ ^ Pi" ■ fl" ■ Tin , 

I'€E„_i /"6E S + 1 

where f € (i?), /// € K[x/'], fp £ -Rp, and //» and rr» are the same as before. 

Notice that all steps in the transformation of the presentation 1)18(1 to the pre- 
sentation (|19J) are invertible. Also if pi ■ fj ■ rj = then all summands in 1)18(1 
and (|19|l coincide. So, if at least the one of the summands in the presentation 1)19(1 
of tfi = is nonzero, then we can also find a nonzero presentation of the form 1)18(1 
for = 0, which is impossible by the induction hypothesis. This shows that the 
presentation 1(19(1 of cj) is unique. 

This proves 1(17(1 . For s — we obtain the claim of Proposition ll0.3l □ 

Finally we can put everything together and prove Theorems 13.11 and 13.31 

Proof of Theorems and Vi.HL For a graph G, let Ag, Bg, and Cg be the algebras 
defined in Sections El El and El Then we have the following termwise inequalities 
of Hilbert series 

(20) Hilb Ac > Hilb B G > Hilb C G ■ 

The first inequality follows from Theorem l5 . 21 because Tq is a monotone monomial 
ideal and, by Lemma 15.11 Jg is its deformation. The second inequality follows 
from Lemma [9.21 that says that Cg is a quotient of Bg- Theorem 12.11 claims that 
dim^G = Nq is the number of spanning trees of the graph G. On the other hand, 
by Theorem l9.ll dimCc = Nq- Thus all inequalities in l(20|) are actually equalities. 
Moreover, by Theorem 19. II the dimensions of fc-th graded components are equal to 

dim^ = dimB G = dimC^ = N l Q l ~ n ~ k , 
the number of spanning trees of G with external activity \G\ — n — k. □ 
Corollary 10.5. The algebras Bg and Cg are isomorphic. 

11. Algebras related to forests 

Definitions of the algebras Bg and Cg and the proof of Theorem 19.11 are similar 
to constructions from PSS1 . Let us briefly review some results from [PSS1 . 

Let G be a graph on the vertices 0, . . . , n. Let Jq be the ideal in . . . , x n ] 

generated by the polynomials 




where I ranges over all nonempty subsets in {1, . . . ,n} and the number Di is the 
same as in SectionEl cf. (@J). Let Bg — K[x±, . . . , x n ]/ Jg- 

Let 4>g be the commutative algebra generated by the variables <p e , e G G, with 
the defining relations: 

((f> e ) 2 — 0, for any edge e. 
And let Cg be the subalgebra of $g generated by the elements 

for i = 1, . . . , n, cf. SectionEl 
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A forest is a graph without cycles. The connected components of a forest are 
trees. A subforest in a graph G is a subgraph F C G without cycles. Fix a linear 
order of edges of G. An edge e £ G \ F is called externally active for a forest 
F if there exists a cycle C in G such that e is the minimal element of C and 
(C \ {e}) C i* 1 . The external activity of -F is the number of externally active edges 
for F. 

Theorem 11.1. |PSSlj The algebras Bq and Cq are isomorphic to each other. 
Their dimension is equal to the number of subforests in the graph G. 

The dimension dimB^j of the k-th graded component of the algebra Bq equals 
the number of subforests F of G with external activity \G\ — \F\ — k. 

In |PSS2| we investigated the algebra Bq for the graph G = K n+ \. Let I n = (rhi) 
and J n = (pi) be two ideals in the polynomial ring K.[x\, . . . , x n ] generated by the 
monomials rhj and the polynomials pi, correspondingly, given by 

pi = (x h + ■ ■ ■ + x lr y^- r +^+\ 

where / — {i\ < - • • < i r }, ranges over nonempty subsets of {1, . . . , n}, cf. Subsec- 
tion Notice that I n is a monotone monomial ideal and j n is its deformation. 
Let A n = K[xi, ... , x n ]/T n and B n = K[xi, . . . , x n ]/J n . 

Let us say that a non-negative integer sequence b = (pi, . . . ,b n ) is an almost 
parking function of size n if the monomial X — X ' ' ' X does not belong to 
the ideal T n . Clearly the class of almost parking functions includes usual parking 
functions. 

For a forest F on the vertices 0, . . . , n, an inversion is a pair of vertices labelled 
i and j such that i > j and the vertex i belong to the path in F that joins the 
vertex j with the minimal vertex in its connected component. 

Theorem 11.2. [PSS11 IPSS2j The algebras A n and B n have the same Hilbert 
series. The dimension of these algebras is equal to the number of forests on n + 1 
vertices. 

Moreover, the dimension dim.A n = dimi?^ of the k-th graded components of the 
algebras A„ and B n is equal to 

(A) the number of almost parking functions b of size n such that Yl^—i bi = k; 

(B) the number of forests on n + 1 vertices with external activity Q) — k; 

(C) the number of forests on n + 1 vertices with (™) — k inversions. 

The images of the monomials x b , where b ranges over almost parking functions 
of size n, form linear bases in both algebras A n and B n . 

Theorem HO first stated in |PSS2j . follows from results of jPSSlj . The algebra 
B n is the algebra generated by curvature forms on the complete flag manifold. It 
was introduced in an attempt to lift Schubert calculus on the level of differential 
forms, see |PSSli rPSS2 ShSh . This example related to Schubert calculus was our 
original motivation. 

12. p- ALGEBRAS AND p-PARKING FUNCTIONS 

We conclude the paper with a discussion of a special class of monotone monomial 
ideals and their deformations. 
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Let p — (pi, . . . , p n ) be a weakly decreasing sequence of non-negative integers, 
called a degree function. Let T p — (mi) and J p — {pi) be the ideals the ring 
. . . , x n ] generated by the monomials m/ and the polynomials pi, correspond- 
ingly, given by 

mi = (x ix ■ ■ -Xi r ) Pr , 

Pi = (a;*! H h Xi r ) r ' Pr , 

where in both cases / = . . . , i r } runs over all nonempty subsets of {1, . . . , n}. 
Let A p = K[x lt . . . ,x n ]/l p and B p = K[xi, . . . ,x n ]/J p . 

Let us say that a non- negative integer sequence b = (b\, . . . , b n ) is a p-parking 
function if the monomial x^ 1 ■ ■ ■ x h ^ does not belong to the ideal I p . More explicitly, 
this condition can be reformulated as follows. A non- negative integer sequence 
b — (b\, ... , b n ) is a p-parking function if and only if, for r = 1, . . . , n, we have 

#{« | b t < Pn-r+l} > r. 

This condition can also be formulated in terms of the increasing rearrangement 
ci < • • ■ < c n of the elements of b as c, < p n +i-i- The p-parking functions were 
studied in JPiSt] and in |Yanj . They also appeared under a different name in |PPj . 
Notice that (n, . . . , l)-parking functions are exactly the usual parking functions of 
size n. 

The monomials x b , where b ranges over p-parking functions, form the standard 
monomial basis of the algebra A p . Thus the Hilbert series of the algebra A p equals 

mhA p = J2<i bl+ '" +K , 

b 

where the sum is over p-parking functions. The dimension dim„4 p of this algebra 
is equal to the number of p-parking functions. 

Theorem 15.21 specializes to the following statement. 

Corollary 12.1. The monomials x b , where b ranges over p-parking functions, lin- 
early span the algebra B p . Thus we have the termwise inequality of Hilbert series: 

mib A p > Hilb B p . 

It would be interesting to describe the class of degree functions p such that 
Hilb^lp = HilbSp. If p r = I + k(n — r) is a linear degree function then, according 
to Corollary 14. 31 the Hilbert series of A p and B p are equal to each other and 

dim A p = dim 6^ = 1(1 + /cn)" -1 . 

For n = 3, Schenck |Schej gave another proof of this fact using ideals of fatpoints. 

Let us say that a degree function p is almost linear if there exists an integer 
k such that pi — Pi+i equals either k or k + 1, for i = 1, . . . , n — 1. Computer 
experiments show that the equality Hilb A p = Hilb B p often holds for almost linear 
degree functions p. The table below lists some almost linear degree functions, for 
which the equality Hilb„4 p = Hilb£> p holds. 

On the other hand, the equality of Hilbert series fails for the almost linear degree 
functions p = (9, 6, 3, 1) and p = (9, 7, 5, 4, 3). We do not know an example when 
Hilb^lp = HilbSp and p is not almost linear. 

The ideal T p is a strictly monotone monomial ideal provided that the degree 
function is strictly decreasing pi > ■ • • > p n > 0. Corollary 16.41 gives a minimal 
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p 


dim„4 p 


P 


dim A p 


(4,2,1) 


25 


(8,5,3) 


306 


(8,5,1) 


142 


(8,6,3) 


351 


(6,4,3) 


153 


(11,7,2) 


506 


(9,5,2) 


290 


(12,8,3) 


855 


(6,4,3,2) 


632 


(8,7,5,3) 


3021 


(9,6,4,2) 


2512 


(9,7,6,5) 


4925 


(8,6,5,3) 


2643 


(11,8,6,3) 


7587 


(8,6,5,4) 


2832 


(12,9,7, 4) 


12460 


(9,8,6,4, 2) 


31472 


(10,9,7,5,3) 


65718 



free resolution for this ideal. Recall that S — K[xi, . . . ,x n ] and S(—d) is the free 
Z-graded S-module of rank 1 with generator of degree d. 

Corollary 12.2. Let p be a degree function such that p\ > ■ • • > p n > 0. The ideal 
X p has a minimal free resolution of the form 

■ ■ ■ — ► C3 — ► C2 — ► C\ — ► Co = S — > S/Xp — > 0, 

with 

c k = s(-d(h,...,i k ))^-- n -'0 , 

h,---,lk 

where the direct sum is over Zi, . . . , Z& > 1 such that h + ■ ■ ■ + Ik <n, 
d(h, ...,l k ) = h p{h) + h p(h + h) + ■ ■ ■ + h p(h + ■■■ + h), 

and (, n , ) = ti — :— 77 — ^ j-tt is the multinomial coefficient. 

\h,— ,l k J hl—l h l(n-h IkY- J 

Coni ecture 16 . 1 Ul imply that if dim^lp = dimS p then the ideals X p and J p have 
the same graded Betti numbers. It is already a nontrivial open problem to prove 
(or disprove) that the graded Betti numbers of the ideal J p , for a linear degree 
function p, are given by the expression in Corollary 112.21 Schenck's results |Sche| 
for n — 3 seem to support this conjecture. 

Proposition l8.1l specializes to an expression for the Hilbert series Hilb A p with al- 
ternating signs. Actually, in this case it is possible to give a simpler subtraction-free 
expression for the Hilbert series. The following statement is a slight enhancement 
of a result of Pitman and Stanley, who gave a formula for the number of p-parking 
functions. 

Proposition 12.3. cf. |PiStl Theorem 11] The Hilbert series of the algebra A„ 
equals 

(2D Hiib^=£n £ — i— - — . 

where the sum is over (n + 1) I1_1 usual parking functions a — (a%, . . . , a n ) of size 
n. Here we assume that p n +i = 0. Thus the dimension of A p , which is the number 
of p-parking functions, is given by the following polynomial in pi, . . . , p n : 

n 

dim A p = ^ Y[(pn- ai - Pn-a, + l), 
a i—1 



where again the sum is over usual parking functions of size n. 
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Proof. For i = 0, . . . , n, let Z% be the interval of integers Zi = [p n -i+i, Pn-i[i where 
we assume that po = +00 and p n +i = 0. Then the set of positive integers is the 
disjoint union of Zq, . . . , Z n . Let / : b t— > a be the map that sends a positive integer 
sequence b = (pi, ... , b n ) to the sequence a = (a±, . . . , a n ) such that 6, € Z Qi for 
i = 1, .. . ,n. Then b is a p-parking function if and only if a is a usual parking 
function of size n. Fix a parking function a of size n. Then 

n 

e ? bi+ - +h " = n e 9*. 

b:f(b)=a i=l bi£Z a . 

is exactly the summand in (|21|l . □ 
For example, the Hilbert series of .A p , for n = 2 and n = 3, are given by 
Hilb -4( Pl , P2 )(g) = [ P 2? + 2^ 2 [pi - p 2 ] [p 2 ] , 

Hilb A( Pl , P2 , P3) (g) - [p 3 ] 3 + 3g^[P2 - p 3 ] [p 3 ? + iq 2p3 [P2 - P3] 2 [ps}+ 

+3q" 2 [pi - p 2 ] [p 3 ] 2 + 6 g " 2 +^ [ Pl - p 2 ] [p 2 - p 3 ] [p 3 ] , 

where [s] = 1 + q + h denotes the ^-analogue of an integer s. 

Finally, we formulate a theorem that gives a combinatorial interpretation of the 
value of the Hilbert series Hilb.4 p at q = —1. This theorem follows from results 
of |PP| on p-parking functions. 

Theorem 12.4. [PP] The number Hilb„4 p (-1) equals the number 

of permutations o~%, . . . , a„ of 1, ... ,n such that 

a x V Pl o~ 2 V P2 •• • V— 1 <t„ V p " 0, 

where the notation a V fc b means that a < b for even k and a > b for odd k. In 
particular, Hilb ^4 P ( — 1) is zero if and only if p n is even. 

This theorem basically says that Hilb _4 P ( — 1) is either or plus/minus the num- 
ber of permutations with prescribed descent positions. 

In the case of usual parking functions of size n, i.e., for p = (n, . . . , 1), this the- 
orem amounts to the well-known result of Kreweras |Krew| that the value of the 
inversion polynomial /„(— 1) = (— 1)W Hilb.4(„ 1) is the number of alter- 
nating permutations of size n. 

13. Appendix: Abelian sandpile model 

In this appendix we discuss the abelian sandpile model, also known as the chip- 
firing game. It was introduced by Dhar |Dhar| and was studied by many au- 
thors. We review the sandpile model for a class of toppling matrices introduced 
by Gabrielov |Gab2| . which is more general than in Dhar . Then we show how 
G-parking functions from Section are related to this model. 

Let A = (Ajj)i<ij<„ be an integer n x n-matrix. We say that A is a toppling 
matrix if it satisfies the following two conditions: 

(22) Aij < 0, for i ^ j; and there exists a vector h > such that A h > 0. 

Here the notation h > means that all coordinates of h are strictly positive. 
Notice that conditions (|22|l imply that An > for any i. These matrices appeared 
in |Gab2j under the name avalanche-finite redistribution matrices. 
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Let us list some properties of toppling matrices. Recall that Lq denotes the 
truncated Laplace matrix that corresponds to a digraph G on the vertices 0, . . . , n, 
see Equation in Section |21 

Proposition 13.1. cf. |Gab2| 1. A matrix A is a toppling matrix if and only if 
its transposed matrix A T is a toppling matrix. 

2. Every integer matrix A such that 

(23) Ay < for i ^ j; A„ > for all i; det A ^ 

3 

is a toppling matrix. Equivalently, the truncated Laplace matrix A = Lq cor- 
responding to a digraph G with at least one oriented spanning tree is a toppling 
matrix. 

3. If A is a toppling matrix then all principal minors of A are strictly positive. 

4. If A is a symmetric integer matrix with non-positive off-diagonal entries, then 
A is a toppling matrix if and only if it is positive-definite. 

Proof. 1. This claim follows from |Gab21 Theorem 1.5]. It also follows from the 
result of |Kacl Theorem 4.3], obtained for classification of generalized Cartan ma- 
trices. 

2. Conditions (|23|l are equivalent to the statement that A = La is the truncated 
Laplace matrix for some digraph G with at least one oriented spanning tree, see 
the Matrix- Tree Theorem, Equation Q in Sectional Let dist(i) be the length of 
the shortest directed path in the digraph G from the vertex i to the root 0, and let 
h(e) = (hi,..., h n ) T , where hi = l- e dist «. Then A h(e) > for sufficiently small 
e > 0. Indeed, the z-th coordinate of the vector A h(e) is 

Oio(l-c)+ ««(e di8t(j ' ) -e di8t(i) ), 
j¥<M 

where the ay are the entries the adjacency matrix of G. The leading term of this 
expression has order of gd^M- 1 and is strictly positive. 

3. The fact that det A > is given in |Gab2l Proposition 1.12]. Let us show that 
it also easily follows from the Matrix- Tree theorem. Let A be a toppling matrix 
and h = (hi, ... , h n ) T > be an integer vector such that A h > 0. Then all row 
sums of the matrix A = A • diag(/ii, . . . , h n ) are positive. This means A = Lq is 
the truncated Laplace matrix for some digraph G. The i-th row sum of A is the 
number of edges in G connecting the vertex i with the root 0. According to the 
Matrix- Tree Theorem, the determinant det A is the number of oriented spanning 
trees in the digraph G. This number is positive, because each vertex is connected 
with the root by an edge in G. Thus det A = (hi ■ ■ ■ hn) -1 det A > 0. Any principal 
minor of A also has positive row sums. The same argument holds for the minors. 

4. This claim follows from |Kacl Lemma 4.5]. □ 

Let us now fix a toppling matrix A, and let A, = (Aji, . . . , Aj„) be the i-th 
row of A. A configuration u — (ui, . . . , u n ) is a vector of non-negative integers. In 
the sandpile model, the number m is interpreted as the the number of particles, or 
grains of sand, at site i = 1, ...,n. A site i is critical if ut > An. A toppling at 
a critical site i consists in subtraction the vector Aj from the vector u. In other 
words, toppling at site i decreases Ui by An particles and increases uj by — Ay 
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particles, for all j ^ i. A configuration u is called stable if no toppling is possible, 
i.e., < m < An for all sites i. 

Dhar |Dhar| assumed that the toppling matrix A has non-negative row sums, 
i.e., he assumed that A = Lq satisfies conditions (|23|l . In this case a toppling 
cannot increase the total number of particles. When a toppling occurs, some of the 
particles at site i are distributed among the neighboring sites and some particles 
are removed from the system. While this condition is important from the physical 
point of view, it is not really necessary for the following algebraic constructions, 
cf. Gabrielov Gab2 . Moreover, there are interesting examples, for which this 
condition fails. We still attribute the following results to Dhar even though we will 
not assume that A has non-negative row sums. The proofs that we include for 
completeness sake are close to the proofs from |Dhar] . 

Lemma 13.2. |Gab2j . cf. |Dhar| Every configuration can be transformed into a 
stable configuration by a sequence of topplings. This stable configuration does not 
depend on the order in which topplings are performed. 

Proof. Conditions 11' I'D for a toppling matrix imply that there exists a vector h > 
such that (h, Aj) > for any i. For a configuration u, the value h(u) is non- 
negative and every toppling strictly decreases this value. Thus, after at most 
h(u)/mm(h, Aj) topplings, the configuration u transforms into a stable configu- 
ration. 

If an unstable configuration u has two critical sites i and j, then j is still a 
critical site for u — Aj. Thus is it possible to perform a toppling at site i followed 
by a toppling at site j producing the configuration u — A, — Aj. This operation 
is symmetric in i and j. Using this argument repeatedly, we deduce that the final 
stable configuration does not depend on the order of topplings. □ 

The avalanche operators A\ , . . . , A n map the set of stable configurations to itself. 
The operator Ai is given by adding 1 particle at site i, i.e., increasing Ui by 1, and 
then performing a sequence of topplings that lead to a new stable configuration. 

Lemma 13.3. |Dhar| The avalanche operators A\, . . . , A n commute pairwise. 

Proof. The stable configuration AiAju is obtained from u by adding a particle at 
site j then performing a sequence of topplings then adding a particle at site i and 
performing another sequence of topplings. If we first add two particles at sites i 
and j, then all topplings in these two sequences are still possible and lead to the 
same stable configuration. This shows that Ai and Aj commute. □ 

The abelian sandpile model is the random walk on the set of stable configurations 
that is given by picking a site i at random with some probability pi > and 
performing the avalanche operator Ai. Informally, we can describe it as the model 
where we drop a grain of sand at a random site and allow the system to settle to a 
stable configuration. 

Dhar described the steady state of this random walk. A stable configuration u 
is called recurrent if there are positive integers Cj such that A^u = u for all i. Let 
7Z denote the set of recurrent configurations. The commutativity of the avalanche 
operators Ai implies that the set 1Z is closed under the action of these operators. 
Moreover, the operators Ai are invertible on the set 1Z. Indeed, A~ 1 u can be defined 
as a1 i ~ 1 u for a recurrent configuration u. According to the theory of Markov chains 
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all recurrent configurations have the same nonzero probability of occurrence in the 
steady state and all non-recurrent configurations have zero probability. 

The sandpile group SG, also known as the critical group, is the finite abelian 
group generated by the avalanche operators Ax,..., A n acting on the set TZ. 

Theorem 13.4. Dhar] The sandpile group is isomorphic to the quotient of the 
integer lattice SG ~ Z™/ (A), where (A) = ZAi © • ■ • © ZA n is the sublattice in 
Z™ spanned by the vectors Aj. The order of this group is equal to the number of 
recurrent configurations and is given by \SG\ = \TZ\ = det A. 

Proof. Since there are finitely many recurrent configurations, we may assume that 
the numbers Ci are the same for all recurrent configurations. For a recurrent config- 
uration u£K and an integer vector v — (vi, . . . ,v n ), let A v u — A^ 1 ■ ■ ■ A^u. Then 
u = A u ~ v v for any u, v £ TZ. Indeed, the configuration A u ~ v v is given by perform- 
ing topplings to the configuration (u — v + Nc) + v = u + Nc, where c = (ci , . . . , c„) 
and N an integer large enough to make the vector u — v + Nc positive. The result 
of these toppling equals A Nc u = u. This shows that SG acts transitively on TZ. If 
an element of SG stabilizes a configuration u £ TZ then, by transitivity, it stabilizes 
any other element of TZ and is the identity in SG. Thus the order of the sandpile 
group SG equals to \TZ\. The bijection between SG and TZ is given by A v i— ► A v - u*, 
where is any fixed element of TZ. 

If we add A^ particles at site i to a configuration u and perform a toppling at 
the (unstable) site i, the result will be the same as adding —Ay particles at all 
other sites j ^ i. Thus 

Af" = Y[ Aj A ' J , or equivalently, A Ai = 1 for any i. 

On the other hand, A v ^ 1 if v (A), since topplings are given by subtraction 
of the vectors Aj and A v u £ u + v + (A). This shows that A" = 1 if and only if 
v £ (A) and the map v i— > A" is an isomorphism between the sandpile group SG 
and the quotient Z™/ZA. 

Finally, the order of IT 1 / (A) equals det A. □ 

Dhar suggested a more explicit characterization of the set TZ of recurrent config- 
urations. Let us say that a configuration u is allowed if for any nonempty subset I 
of sites there exists j £ I such that 

Proposition 13.5. |Dharj Every recurrent configuration is allowed. 

Proof. Let u be a recurrent configuration. Then A c u = u, where c = (ci, . . . , c„) 
and Ci > for all i. This means that there exists a sequence of sits i\, . . . ,ik such 
that (i) Ajj + • ■ ■ + A ik = c; and (ii) u + A^ + • • • + A ir > for any r = 1, . . . , k. 
Since all coordinates of Aj, except the i-ih coordinate, are non-positive and c > 0, 
condition (i) implies that the sequence ix, ■ ■ ■ ,ik contains all sites 1, . . . , n at least 
once. 

Let us say that a configuration v is I-forbidden, for some subset / of sites, if 

0<Vj< ]T (-Ay), 
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for all j £ I. If v is /-forbidden and v+Ai > 0, for some site i, then the configuration 
v + Ai is J\ {i}-forbidden. Also notice that there are no 0-forbidden configurations. 

Suppose that the recurrent configuration u is not allowed. Then u is /-forbidden 
of some subset /. We obtain by induction on r = 0, . . . , k that the configuration 
it(, r ) — u + Ajj + • • • + Ai r is / r -forbidden, where I r = I \ {i\, . . . , i r }. In partic- 
ular, it( fe ) is //.-forbidden, where Ik — 0, which is impossible. This shows that the 
configuration u is allowed. □ 

Dhar suggested that a configuration is recurrent if and only if it is stable and 
allowed. Gabrielov |Gabll Section 3, Appendix E] showed that this statement is 
not true in general, and proved the conjecture for a toppling matrix A with non- 
negative column sums and, in particular, for a symmetric toppling matrix A = Lq 
corresponding to an undirected graph G. For symmetric A — Lq, Dhar's conjecture 
was also proved by Ivashkevich and Priezzhev IvPr], and recently by Meester, 
Redig, Znamenski MRZ] Theorem 5.4], and by Cori, Rossin, and Salvy |CRSI 
Theorem 15]. 

The following two claims show how G-parking functions from Section [21 are re- 
lated to the sandpile model. For a vector u = (u\, . . . , u n ), let u v = {u{, . . . , u^), 
where = An — 1 — Uj. 

Lemma 13.6. Let G be a digraph with at least one oriented spanning tree, and let 
A = Lq be the transpose of the truncated Laplace matrix for the digraph G. For 
the sandpile model associated with the toppling matrix A, a configuration u is stable 
and allowed if and only if u v is a G-parking function. 

Proof. Parts 1 and 2 of Proposition 113.11 imply that A = Lq is a toppling ma- 
trix. The statement of the lemma is immediate from the definitions of allowed 
configurations and G-parking functions. □ 

A toppling matrix A is the transpose Lq of the truncated Laplace matrix for 
some digraph G if and only if it has non-negative column sums: 

^2 &ij > for an y 3- 

i 

Theorem 12 . II recovers Gabrielov's result on recurrent configurations. 

Corollary 13.7. [Gabll Eq. (21)] For a toppling matrix A with non-negative 
column sums, a configuration is recurrent if and only if it is stable and allowed. 

Equivalently, a configuration u is recurrent if and only if u v is a G-parking 
function, for G and A — Lq such as in Lemma \V3.b\ 

Proof. Theorem 12.11 and the Matrix- Tree Theorem imply that the number of stable 
allowed configurations equals to dct Lq = dot A. According to Theorem 113.41 the 
number of recurrent configurations is also equal to det A. These facts, together 
with Proposition 113. 51 imply the statement. □ 

Remark that we need to impose the transpose of Dhar's physical conditions 123|) 
on the toppling matrix A in Corollarv ll3.7l The number of recurrent configurations 
for a toppling matrix A is equal to the number of recurrent configurations for the 
transposed toppling matrix A T , because det A = det A T . It would be interesting 
to present an explicit bijection between these two sets of configurations. 
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